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ABSTRACT 

The mechanism responsible for the natal kicks of neutron stars continues to be a challenging 
^\J I problem. Indeed, many mechanisms have been suggested, and one hydrodynamic mechanism 

may require large initial asymmetries in the cores of supernova progenitor stars. Goldreich, 
Lai, & Sahrling (1997) suggested that unstable g-modes trapped in the iron (Fe) core by the 
convective burning layers and excited by the e-mechanism may provide the requisite asymmetries. 
We perform a modal analysis of the last minutes before collapse of published core structures 
{/^ , and derive eigenfrequencies and eigenfunctions, including the nonadiabatic effects of growth by 

^O ' nuclear burning and decay by both neutrino and acoustic losses. In general, we find two types 

^"T ' of g-modes: inner-core g-modes, which are stabilized by neutrino losses and outer-core g-modes 

f^ , which are trapped near the burning shells and can be unstable. Without exception, we find at 

least one unstable g-mode for each progenitor in the entire mass range we consider, 11 M© to 
O ,■ 40 Mq. More importantly, we find that the timescales for growth and decay are an order of 

magnitude or more longer than the time until the commencement of core collapse. We conclude 
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Vh ' that the e-mechanism may not have enough time to significantly amplify core g-modes prior to 

collapse. 
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1. Introduction 

Observations indicate that many pulsars have high proper motions (Lyne & Lorimer 1994; Lorimer, 
Bailes, & Harrison 1997; Hansen & Phinney 1997; Cordes & Chernoff 1998) and that many neutron star 
systems require kicks at birth. For example, pulsar bow shock morphologies (Cordes, Romani, & Lundgren 
1993), misaligned spins of binary pulsars (Cordes, Wasserman, & Blaskiewicz 1990; Kramer 1998; Wex, 
Kalogera, & Kramer 2000; Kaspi et al. 1996; Lai, Bildsten, & Kaspi 1995; Lai 1996; Kumar & Quataert 
1997), the statistics of low-mass X-ray binaries (LMXBs) (Kalogera 1997; Kalogera & Webbink 1998), the 
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high orbital eccentricities of some Be star systems (Verbunt & van den Heuvel 1995), and scenarios for double 
neutron star systems (Dewey & Cordes 1987; Fryer & Kalogcra 1997; Fryer, Burrows, & Benz 1998; Dewi 
& van den Heuvel 2004; Willems & Kalogera 2004) all indicate that neutron stars experience significant 
natal kicks. Furthermore, Arzoumanian. Chernoff. & Cordes (2002) find a bimodal distribution of birth 
kick speeds with peaks at ~100 km s^^ and ~700 km s^^ and dispersions of 90 and 500 km s^^. They 
suggest that roughly 50% of the pulsars in the solar neigborhood escape the Galaxy and that ~15% have 
kick speeds greater than 1000 km s~^. The observations indicate that kicks seem to be hallmarks of neutron 
star formation. 

What is the mechanism producing the observed and inferred kicks? To give a 1.4 Mq neutron star 
a speed of 1000 km s~^, the mechanism must be able to impart ~10^^ erg of kinetic energy, which is 
just 1% of the canonical SN mechanical energy (~10^^ erg) or 0.01% of the collapse energy (~10^^ erg). 
Additionally, any explanation must address the distribution of speeds. If the distribution is bimodal, then 
it is possible that more than one mechanism is responsible for the kicks. Lai (2000b) has summarized the 
possible magnetic, neutrino (i'), and hydrodynamic kick mechanisms, but in general many of the scenarios 
have trouble producing kicks larger than ~200 km s~^. 

Of the many mechanisms propounded, hydrodynamic mechanisms are the most compelling. For exam- 
ple, neutrino-driven convection, which is inherently a nonspherical phenomenon, excites large-scale convective 
motions that stochastically jostle the protoneutron star (Burrows, Hayes, & Fryxell 1995; Janka & Miiller 
1994). However, multi-dimensional calculations indicate that any kick imparted as a result of the convective 
jostling ("Brownian motion") has an average magnitude of ^-^200 km s~^ (Burrows et al. 1995; Janka & 
Miiller 1994; Scheck et al. 2004) and cannot explain the highest velocities nor the second peak in a possible 
bimodal distribution. 

These analyses suggest the need for additional asymmetries before the onset of collapse. Burrows & 
Hayes (1996) state that global asymmetries in the progenitor may be prerequisites for large kicks. In a 
2-D hydrodynamic simulation and with modest asymmetries. Burrows & Hayes (1996) are able to produce 
a recoil speed of ~550 km s^^. However, these calculations are crude and the results, while suggestive, 
are not conclusive. Goldreich et al. (1997) have proposed that these asymmetries may be excited by the 
e- mechanism in which g-mode oscillations in the core are pumped by nuclear reactions. 

As a prelude to more definitive multi-dimensional radiation hydrodynamic investigations of neutron star 
kicks in the supernova context, one would like to determine whether global asymmetries prior to core collapse 
in fact occur. Hence, a thorough analysis of the progenitor oscillations and their growth rates is necessary. 
In this paper, we explore the viability of the g-mode oscillation origin for putative pre-coUapse anisotropics 
by performing such an eigenmode analysis. 



2. Linear Perturbation Equations 

To set the stage for our eigenmode analysis of the cores of nonrotating massive stars, we first describe 
linear stellar pulsation theory in the general case. In linear perturbation theory, a star is assumed to be 
in hydrostatic equilibrium, around which perturbations are small, making second-order or higher terms 
in the perturbation equations negligibly small. Therefore, if fo{r) is a quantity satisfying the hydrostatic 
equations, then f{r^t) -^ fo{f) + f'ir^t), where denotes the background solution and ' denotes the Eulerian 
perturbation. We assume that the perturbation is of the form /'(r, t) = f'{r)Yim{0, ^)e*'^*, where Yi,n{(^, (/>) is 
a spherical harmonic in that i and m are the polar and azimuthal indeccs, repsectively, and lu is the oscillation 



frequency. Under these assumptions and an assumption of adiabatic perturbations, if P, p, 0, Cg are the 
pressure, density, gravitational potential, and sound speed, respectively, the adiabatic linear peturbation 
equations^ are 

i^ + 4^' + (iV^--% = -f^, (2) 

p dr pcj. ar 

r^ dr \ dr J r^ VPC^ g 



and 



r'^ ~ -^-T^^' - 47rGp _ + _e. . (3) 



^r and ^h are the radial and horizontal components of the Lagrangian displacement vector, ^ = r — ro, and 
are defined by 

^^= U,^H§-^,^n-^) YUOA)e'-' (4) 



and 

a = ^f- + 0'V (5) 

^ ' \ P J 

N and L^ are the Brunt- Vaisala and Lamb frequencies, respectively, which are given by 



and 

^t — 2 ' 

where g = GMr/r"^ , G is the gravitational constant, and Mr is the mass interior to r. 



(7) 



Equations (1), (2), and (3) represent a fourth-order boundary value problem. Near r ~ 0, simplifications 
and regularity in the variables (/>', S,r, and [P' / p + 4>') dictate that the inner boundary conditions are 

^-^ = (8) 

dr r 

and 

^r - e^h = 0. (9) 

At the surface of the calculational domain, r = R or Mr = M, one outer boundary condition, from the 
regularity of 0' , is 

^ + (l±ll0' = 0. (10) 

dr r 

A second outer boundary condition allows for the possibility that the atmosphere is either evanescent to 

both p- and g- waves or progressive to outgoing p- waves (Lai 2000b), depending upon the mode frequency: 

'A_-(T//ri-3)\ ^r u;\ /ai(A_-(y/ri-3)) 



Oi J r g \ 6i J gr 

where 



^For a complete derivation of the adiabatic linear peturbation equations sec Unno et al. (1989). 



9 


^' GM Ti' 

'^ - GM g ' 

_ 6i(r7VV5) - ((1 + rN^g) + i)V/ri 


9 



(13) 
(14) 
(15) 

"'^ {V/ri-3 + i){l + rNyg + e)-hb2' ^^^^ 

and 

^ b2V/T,-iV/T,-3 + e)rNyg 

' iV/ri-3 + l)il + rNyg + e)^bib2' ^ ' 

If the mode frequency is such that this boundary condition represents leakage of mode energy by progressive 
p- waves then the eigcnsolutions arc complex, with the imaginary part of w, w/, being related to the rate of 
leakage. 

The solutions to the boundary and eigenvalue problem as defined by eqs. (l)-(3) with the boundary 
conditions, eqs. (8)-(ll), are the eigenmodes of the perturbed star. The eigenvalue, w^, and eigensolution 
depend on the stellar structure and spherical harmonic index, €, but without rotation are degenerate for 
different values of m. We have found solutions to this boundary value problem using a shooting technique, 
in which we search for all modes within a region of the complex plane of w^. This is accomplished with a 2-D 
Newton- Raphson technique, where convergence to one part in 10^, or better, in the eigenvalue is required. 
We compare our eigensolutions and eigenfrequencies for the standard Solar model (Christensen-Dalsgaard 
et al. 1996) with the eigenmodes determined by a publicly accessible Fortran code^ designed to solve the 
adiabatic linear perturbation equations. The agreement between the two is better than 0.1%. 



2.1. Work Integral 

In general, non-adiabatic effects such as nuclear burning and neutrino emission are included in the 
perturbation equations. However, if these terms are small, the contributions to the imaginary part of w, w/, 
which is directly related to the rate of growth or decay of the mode, can be obtained from the time-averaged 
work integral, Wj (Unno ct al. 1989). For a quasi-adiabatic oscillation. 






■' c.7o " V"^ r3-i 



^-^YUeA) 



dMr, (18) 



where ej (erg g^^ s^^) is the deposition or loss of energy by either nuclear heating if j = nuc or neutrino 
losses if j = v. Furthermore, ejx ~ {dlnej/dlnT)p, Sjp — {dlnej / dlnp)T, and Fa — 1 = {dlnT/dp)s, where 
S is the entropy. The radial distribution of the Lagrangian perturbation in temperature is given by 

^.v'v..(^a-^-^), (19, 

T \ g gr r J 

where Vad ~ [dlnT /dlnp)s and V = grp/P. 

^provided by J. Christensen-Dalsgaard http://astro.phys. au.dk/~jcd/adipack.n/ 



When Lo is complex, as is the case when the mode gains or loses energy, the time dependence of the 
modes may be written as 

^"oce'"'«*e*/^, (20) 

where uju is the real part of w, and r is the timescalc for growth (+) or decay (-) and is a composite effect 
determined by the equation 

1.^ + 1 + ^^ (21) 

Tnuc, Ty, and Tipak includc the effects of nuclear gains, neutrino losses, and acoustic losses, respectively. The 
nuclear and neutrino timescales are related to their respective work integrals by 






where Ew is the energy in the mode: 



p — fl 



c? '■^' 



dMr, (23) 
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and Tieak is directly determined by the imaginary part of lu, Im{ui) = —l/rieak- This quantity is the solution 
to the adiabatic linear perturbation equations when eq. (11) is the fourth boundary condition. The modes 
are stable or unstable for r < or r > 0, respectively. 



3. Massive Star Structures 

To investigate the character of the eigenmodes as a function of progenitor mass and time, i, prior to 
the onset of collapse, we consider the cores of nonrotating, massive stars with masses of 11 M©, 13 M©, 15 
Mq, 25 Mq, 30 Mq, 35 M© and 40 M© (Rauscher et al. 2002; Woosley, Heger, & Weaver 2002) and at 
t = 3600, 600, 60, 30, 25, 20, 15, 10, 5, 2, 1, and (start of collapse) seconds before the onset of collapse. 
The commencement of collapse is defined at the time when the peak infall velocity exceeds 900 km s~^. 
This is about 200 to 250 ms before core bounce. The structure of each stellar model uniquely determines 
its eigenmodes. Thus, an indcntification of important trends in the models is helpful in characterizing these 
modes. In the top panel of Fig. 1 we see that more massive models have higher entropy, translating into lower 
central densities and more extended density (p) versus interior mass, Mr, profiles. Consequently, prominent 
boundaries which are associated with discontinuities in entropy and Ye and which affect the character of the 
modes are also located at larger radii for higher-mass progenitors. These discontinuities are the boundary 
between the fossil Fe core (the residue of core Si burning) and the most recent Fe ashes from shell Si burning 
(circles on Fig. 1), the edge of the Fe core (includes the fossil Fe core and recent Fe ashes) and the Si burning 
shell (stars on Fig. 1), and the interface between the Si burning and O burning shells (triangles on Fig. 1). 
The general structural profiles and these distinct boundaries (Table 1) together determine the eigenfunctions 
and eigenfrequencies. 

The role of structure in affecting the general character of the solutions can be understood using the 
local dispersion relation for plane waves of the form e'(*''''^"*): 



'"'-"^'f-''. (24) 



where fc^ is the radial component of the wave vector (Unno et al. 1989). This relation implies that for 

-2 



N'^ < u!^ < Lj, all waves are evanescent, for N'^ < Lj < uP' , they are p- waves, and for lo^ < N"^ < Lj, 



they arc g-wavcs (Unno ct al. 1989). A propagation diagram, which plots N'^ and L^ versus Mr or r, shows 
where in the star each type of wave is allowed to propagate. The propagation diagrams for the 13 Mq (top 
panel) and 30 M© (bottom panel) models and their corresponding density and Y,, profiles, are presented 
in Fig. 2. The positions of the prominent discontinuities are marked as in Fig. 1. Rather than plotting 
N"^ and Lf, the equivalent period has been presented for each. These periods are Psrunt = '2t'/^/\N^ and 
Phamb = 27r/-\/Lj. An important feature illustrated in Fig. 2 is that A''^ rises and then falls, creating 
a resonant cavity for g- modes. The fall is partly due to the density and pressure structures, but on the 
outside it dramatically plummets due to convection in the 0-burning shell, thereby confining g-waves within 
the core. For higher frequencies (lower periods), p- waves are allowed. However, they are free to propagate 
throughout the entire star. Since the sound-crossing time in the outer stellar mantle is extremely long and 
dissipative processes exist, the establishment of standing p-modes is prevented. Hence, for all models and 
times, g-modes are trapped in the core, encompassing the Fc core and frequently the Si burning shell. 

Other salient features evident in Fig. 2 are the spikes in PBrunt, or equivalently, A^^. These spikes are 
the result of S, Y^, and p discontinuities at the boundaries indicated in Table 1. As impedance mismatches 
these spikes effectively reflect g-waves. Therefore, even though most g-modes have their greatest amplitude 
near the center of the core, the reflecting boundaries imply that there might be sub-regions within the overall 
resonating cavity in which modes may be trapped. 

By prodigious neutrino emission, the cores of massive stars continue to change their structures prior to 
the onset of core collapse (Woosley et al. 2002). Figure 3 illustrates the changes of the propagation diagram 
for the 13 M© model during the last ten minutes of evolution. As the inner regions quickly evolve and 
contract, the Brunt- Vaisala frequency increases with time, gradually increasing the modal frequencies. In 
contrast, the model corresponding to f = s (onset of collapse) exhibits a dramatic truncation in the size of 
the resonating cavity, in that N"^ drops to negative values as the convective region grows due to implosive 
burning. However, collapse happens on a dynamical timescale, and will take 200 to 250 ms. Given that the 
convective turnover times and the periods of oscillation are longer than this, the assumption of steady-state 
mixing-length convection is probably violated. Although for completeness, results for the f = s models are 
presented here, the quantitative results at t = s should be viewed with healthy skepticism. At other times, 
the structural changes are gradual enough for our peturbation theory to be reliable. 

In this paper, the major sources and sinks of energy considered are nuclear burning, neutrino losses, and 
losses by coupling to outward propagating p- waves (Lai 2000b). Since losses by radiation is overwhelmed 
by neutrino emission, perturbations in the radiative flux is ignored. Throughout the Fe core, neutrino losses 
dominate over energy deposition. Only in the Si- and 0-burning shells do nuclear reactions allow the e- 
mechanism to operate. For the typical g-mode with its largest amplitudes in the center, neutrino losses will 
dominate in eq. (21), resulting in stability. However, potentially unstable modes are those that are trapped 
mostly in the outer part of the core near or in the Si-burning region by the discontinuities in entropy. Ye, 
and density (Lai 2000b). 



4. Results 

The results of our analysis are summarized in Table 2, which lists all unstable modes found, and Figs. 
4-11, which are described in detail below. Even though we did find modes for £ = 1, 2, 3, and higher £s, 
we present the results only for £ = 1, 2, and 3 in Table 2. since the trends at higher £s are adequatly 
demonstrated with the set presented in this paper. In general, the most interesting spherical harmonics to 



consider arc those with odd £, or those with a top-bottom asymmetry. However even £s are interesting as 
well, as the nonlinear evolution of the perturbations into a kick may be triggered by any low order L Since 
the dipole, £ = 1, spherical harmonic most closely resembles the desired kick geometry, we focus in this 
section on ^ = 1. with peripheral mention of ^ = 2 and 3. 

Stable and unstable modes are found in all progenitor models. Figures 4 and 5 display the propagation 
diagram for the 13 Mq and 30 Mq models, at t = 5 s before the commencement of collapse. These diagrams 
represent modes of the compact and extended profiles, respectively. Note again, the spikes resulting from 
Yf, and p discontinuities in PBrunt and the existence of a resonance cavity confined by O-shcU convection 
in the 13 Mq model (Fig. 4). Because the O shell is farther out in mass for the more massive 30 M© 
model (Fig. 5), the density and pressure profiles are more relevant in defining the resonant region. In 
both figures each horizontal line corresponds to a specific mode; the shading encodes |^p. Dark and lighter 
shading signify low and high displacement, respectively, while blues denote stable modes and reds identify 
the unstable modes. Figures 4 and 5 show that there are many more stable modes than unstable modes, 
which is true for all models at all times and for all Is. We also identify two classes of g-modes, inner-core and 
outer-core g-modes. Inner-core g-modes have their greatest displacements toward the center and diminish 
with radius. Their spacing in log j^g (Period) systematically decreases with increasing period and increasing 
number of half- wavelengths, n. In contrast outer-core g-modes do not have their greatest amplitude in the 
center. They are trapped by discontinuites between the fossil Fe core and the O shell in the case of the 13 
Mq model and within the Si burning shell in the case of the 30 Mq model. Similar trapping exists for all 
other progenitors we consider. 

Less than half of the outer-core g-modes are modes trapped entirely under the spikes (i.e. surface g- 
modes). While these modes satisfy the adiabatic linear oscillation equations, their physical significance is not 
well constrained. Their eigenfunctions and eigenfrequencies depend on the spike widths and heights, which in 
turn depend upon the interface between convective and radiative regions in mixing-length theory. However, 
the results and conclusions we present in this paper are not predicated on their character or existence. 

A sample of inner-core g-modes for the 13 Mq, t = 5 s, model with ^ = 1 is plotted versus interior mass 
in Fig. 6. The decreasing envelope of the eigenfunction and systematic increase in period with increasing 
half- wavelengths, n, is well explained by a WKB analysis (Unno et al. 1989). The quantization condition, 
WK = J kj-dr, integrated over an appropriate resonant cavity, for g-waves gives 

^,jm2i!if ^*. (25) 

'^^ J cavity ^ 

Numerically, the inner core frequencies are proportional to [£{£ + l)]^"/ri,, and integrating eq. (25) over 
the appropriate resonance cavity gives approximately the same value. Following the conservation of wave 
fiux arguments in Goldreich & Wu (1999), the envelope for the aribtrarily scaled amplitudes of £,r/R and 
^h/R in Fig. 6 also correspond with those from the WKB analysis. In contrast, a simple WKB analysis 
fails to predict either the shape or the period of the outer-core g-mode shown in Fig 7. Instead, this mode's 
period and large displacement farther out in mass are determined by both resonance in the overall cavity 
and resonant trapping between disontinuities in density. Ye , and entropy. 

In the context of losses and gains in energy. Fig. 8 displays |^p versus interior mass for two representative 
inner-core g-modes (top panel) and an outer-core g-mode (bottom panel) for the 13 Mq, t = 5 s model for 
^ = 1, in addition to enuc and e,y in erg g~^ s~^. Figure 9 presents similar information for the 30 Mq model. 
The overlap of [^p with e„„c and e^ is an excellent diagnostic proxy for the integrand in the work integral 
(eq. (18)). Clearly, neutrino losses dominate the work integral for the case of inner-core g-modes in both 



models. For the particular outer-core g-modes plotted for the 13 M0 and 30 A/q progenitors in the bottom 
panels of Figs. 8 and 9, respectively, the large amplitude in the Si-burning region allows for instability. 
While all outer-core g-modes in Fig. 5 are unstable, two of the three 13 Mq outer-core g-modes in Fig. 4 
are stablized by neutrino and acoustic losses. In general, at least one outer-core g-mode has the requisite 
radial perturbation distribution to be unstable. 

The growth timescales versus period for ^=1,2, and 3, t ~ 60, 30, 25, 20, 15, 10, 5, 2. 1, and s before 
the onset of collapse are plotted for 11 Mq 13 Mq, and 30 Mq models in Figs. 12, 10, and 11, respectively. 
As expected, these figures illustrate that for higher £ the periods are smaller. Even though we find many 
unstable modes for the 11 Mq we do not find clear trends in the r versus period plane. In contrast, as 
represented by the 13 Mq and 30 Mq models, all other masses we consider demonstrate a distinct trend. 
The arrows in Figs. 10 and 11 indicate the direction in which the modes evolve with time in the r versus 
period plane. Specifically, later models are more compact and have shorter period oscillations. In addition, 
the later, more compact, models have slightly higher densities and temperatures in the burning regions, 
thereby dramatically increasing the burning rates and decreasing r in Figs. 10 and 11. Also note that at 
no point in any of the models are the growth timescales shorter than the times until the start of collapse. 
These trends continue to the onset of collapse, which is marked on Figs. 10 and 11 with red symbols. In 
fact, for the latest models, t = s and 1 s before the start of collapse, r becomes comparable to the period, 
calling into question the assumption in §2.1 of small nonadiabatic effects. While dramatic changes imposed 
by dynamical collapse and implosive burning for the i := s models make the effects listed above more 
pronounced, the qualitative trends highlighted still obtain for the t = s models. 

For all unstable modes found Table 2 lists the progenitor mass, time until onset of collapse, the period, 
growth timescale, r, and the contributions, t„uc, t^-, and Tieak- From tens of minutes before commencement 
of collapse until collapse itself, unstable modes exist in all progenitor models. However, at every evolutionary 
stage the timescale for growth is an order-of-magnitude or two longer than the time until onset of collapse. 
The rates of growth or decay we derive are too long compared to the time to onset of collapse to significantly 
alter the mode amplitudes. It seems that even though the e-mechanism is present in the later evolutionary 
stages of the cores of massive stars, the timescales imply that its effect is not large. 



4.1. Analytic Estimates 

Consistently, we find that the timescales for growth are longer than the time until the onset of collapse. 
As such the e-mechanism is not able to significantly affect the amplitudes of the unstable g-modes. Hence, our 
evaluation of the succcsfulness of the e-mechanism depends upon the validity of the timescale calculations. 
In this section we corroborate, with an analytic estimate for Tnuc, our numerical evaluations for the growth 
timescales. 

Wj, given by eq. (18), is the rate of energy gained or lost, and Ew, cq. (23), is the energy in the mode. 
Assuming the region where nuclear burning interacts with a growing mode is small and the mode itself may 
be represented by a top-hat function located entirely and only within the burning region with maximum 
displacement Ar/r, then the integrand in eq. (18) is roughly constant. Therefore, 

W, ^ -e,e,. i-\ ^, (26) 



and cq. (19), simplifies to 

f->'V,.(^)^ ,27, 

As a harmonic oscillation, g- modes are the oscillations of isodensity contours about their equilibrium positions 
in hydrostatic equilibrium. The harmonic potential is t/ ~ i9^0e///9r^a;^, where x = Are'"*. From 
V(/) = SIPj p, we get (^f^fj ^ (f) r^ P/ p ~ GMr/r. Hence, the potential energy, u, is given by 



GM, 

u ^ 



-(7) ■ ('') 



r \r - 
Since Ew = 2 J < u> dMr, where <> denotes the time average over one period of oscillation, then 

e^.2Ml(^)'^^ ,2,, 

r \ r J Att 
Finally, taking the ratio of eq. (29) and eq. (26) gives a rough estimate for r: 



'JttEw ^ P/jpV) 



Tnuc = — f77 4- -5-, (30) 



In the burning region, for the 30 M© t = 5 s model and the 8.67 s, ^ = 1 mode, P/ p ~ 2.2 x 10^^ erg g~^, 
f-nucf-nucT ~ 2 X 10^^ erg g^^ s^^, Vad ~ 1/4, and V ^ 2, giving Tnuc ^ 320 s from eq. (30). The numerical 
solution in Table 2 is Tnuc = 328 s. Similarly for the 8.45 s, £ = 1 mode of the 13 M© progenitor model at 
i = 1 s, P/p ~ 1.4 X lO^'^ erg g"\ Cnuc^nucT ^ 10^^ erg g"^ s"\ V ad ~ 1/4, and F ~ 3, giving t„„c ~ 30 
s from eq.(30). The numerical solution (Table 2) is Tnuc = 15.5 s. That the simple analytic arguments 
and estimates approximately reproduce the numerical results is reassuring. Even though the nuclear rates 
increase toward later times, at no time is this trend enough to cause significant growth prior to collapse. 



5. Summary and Conclusions 

For all progenitors, we find stable and unstable g-modes with oscillation periods in the range from ^^1 
to ~10 seconds. The stable g-modes are often concentrated in the inner core and are stabilized by neutrino 
emission. Unstable g-modes on the other hand are trapped in the outer radii of the core by discontinuities 
in 5*, density, and Y^. Their typical growth timescales determined numerically and analytically range from 
10s to 10,000s of seconds, which arc long compared to the time until the start of collapse. Therefore, we 
conclude that the e-mechanism is an unlikely source for large perturbations in the progenitor prior to the 
onset of coUapse. 

There are several caveats to our conclusion. Our results are for a set of 1-D nonrotating models from 
Rauscher et al. (2002) and Woosley et al. (2002). Rotation in the late stages of massive stellar evolution may 
significantly effect the structure of the star prior to collapse (Heger, Langer, & Woosley 2000; Heger, Woosley, 
& Langer 2003), thereby altering the analysis in this paper. However, with modest Fc core rotation periods 
(> 10 s) we expect little deviation from our results, and the simple energy arguments in §4.1 should hold. 
We also have not considered a different suite of progenitor models (Nomoto & Hashimoto 1988; Limongi & 
Chieffi 2003), which may have different eigenmodes. In addition, calculating full 3-D, non-linear, dynamical 
convection with nuclear reactions during the final stages of massive star evolution is a task for the future; 
the progenitors and their g-modes may be significantly different than assumed or calculated here. 
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Wc also have not considered coupling between convectivc modes and g-modes, which provides an excel- 
lent mechanism for the excitation of the 5-minute p- modes in the Sun (Goldreich & Kumar 1990). In such 
analyses, the amplitudes of the modes are estimated by assuming steady state is achieved in the nonlinear 
time-dependent amplitude equations. However, the short timescales until the onset of collapse would violate 
the assumption of steady state and invalidate such analyses. More importantly, it would be unlikely that the 
short timescales until the onset of collapse would allow for significant excitation of g-modcs by convection. 

On the other hand, the growth of perturbations in the supersonic regions during collapse does occur 
(Lai 2000a; Lai & Goldreich 2000). Because we have not ruled out all possible sources of perturbations 
prior to collapse, asymmetries in the progenitor may still be relevant in producing the largest kicks. For 
example, convection in the last stages of evolution is often quite vigorous and may have Mach numbers in 
the range 0.1 to 0.2 and density perturbations, 5p/ p, in the range 0.01 to 0.05 (Bazan & Arnett 1998). These 
perturbations may indeed be sufRcient to produce the observed kicks (Burrows & Hayes 1996). 

In performing a stellar pulsational and stability analysis of the late-time cores of massive stars, we 
have identified many stable and unstable modes. We note that stability is predominantly determined by 
the neutrino losses or energy deposition by nuclear reactions. Compared with the time until the onset of 
collapse, the timescales for growth and decay are long. In §4.1, we argue that the long growth timescales 
(compared with the time until the start of collapse) are confirmed by energetic arguments. We suggest that 
the e-mechanism for g-mode growth does not generate the requisite perturbations prior to collapse needed 
to stimulate large kicks during the supernova explosion. Therefore, we conclude that if the hydrodynamic 
mechanism in combination with progenitor perturbations is to succeed in explaining neutron star kicks other 
sources of seed perturbations must be investigated. 
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this work under the auspices of the U.S. Department of Energy at the Los Alamos National Laboratory 
operated by the University of California under contract No. W-7405-ENG-36. The LANL report number 
for this article is LA-UR-04-3513. J.W.M. would like to thank the Joint Institute for Nuclear Astrophysics 
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Table 1. Structural Boundaries 



Mass (Mq)1 


Fossil Core (Mq)^ 


Fe Core (Mo)^ 


Base (Mq)^ 


11 


0.83 


1.22 (C) 


1.36(C) 


13 


0.98 


1.33 (R) 


1.44(C) 


15 


1.05 


1.41 (R) 


1.74(C) 


20 


0.99 


1.46 (R) 


1.60(C) 


30 


1.06 


1.45 (R) 


2.03(C) 


35 


0.96 


1.47 (M) 


1.66(C) 


40 


0.97 


1.54 (R) 


1.76(C) 



The initial mass of the model. 

^The location of the fossil Fc core, which was formed by convective 
core Si burning, at 5 s before collapse. 

^The location of the base of the Si burning shell or the top of the 
Fe core, which inludes the fossil Fe core and the more recent ashes 
from the Si burning shell. A C indicates that the shell is convective, 
R indicates that it is radiative, and M inidicates that the shell is 
radiative in its inner portion and convective in its outer portion. 

^The location of the base of the O burning shell. 
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Table 2. Growing Modes 



Mass (Mo)i t (s)2 Period (s)^ t (s)"* 1/t„„c (s-1)'5 -1/t, (s-i)*^ -1/rieafe (s-^)^ 



11 2.21 1.23 X IQi 1.60 x 10"^ 6.66 x lO'^ 1.23 x 10"^ 

1 7.63 5.93 x 10^ 6.62 x lO'^ 4.93 x 10"^ ^ o 

1 4.37 1.20 X 10^ 1.52 x 10"^ 6.87 x lO'^ 2.20 x lO^^ 

1 3.04 5.49 x 10^ 2.87 x 10"^ 6.40 x 10"^ 4.12 x lO^^ 

1 2.74 1.86 x 10^ 2.56 x 10"^ 6.30 x 10"^ 1.39 x 10^^ 

5 4.97 1.85 X 10^ 7.09 x lO^^ 1.68 x 10"^ 1.56 x 10"^ 

5 3.53 2.56 x 10^ 7.96 x lO'^ 2.02 x 10"^ 2.03 x lO'^ 

30 9.14 3.47 x 10^ 3.79 x lO'^ 9.07 x 10"^ ~ 

30 6.24 5.69 x 10^ 3.07x10^3 1.31 x lO'^ ~0 

30 5.75 9.77 x 10^ 1.11 x 10"^ 8.26 x lO"" - 

30 4.96 8.20 x 10^ 2.68 x lO'^ 1.44 x lO'^ 1.87 x lO'^ 

30 4.09 1.83 x 10^ 7.23 x 10"^ 9.63 x lO"'' 8.07 x 10"^ 

30 3.62 1.30 x 10^ 4.14 x lO'^ 1.15 x lO'^ 2.22 x 10-^ 

60 9.80 6.74 x 10^ 2.09 x lO'^ 6.11 x lO"'' ~ 

60 6.80 7.24 x 10^ 2.19 x 10^^ ^n ^ 19-4 ^ q 

60 6.37 2.34 x 10^ 4.89 x 10"^ 6.15 x lO"'' ~ 

60 4.48 8.15 X 10^ 2.36 x 10"^ 6.85 x lO"'' 4.50 x lO""* 

13 4.44 8.59 x 10° 2.54 x lO'^ 1.37 x lO'^ 6.46 x 10"^ 

4.34 4.70 X 10-1 2.15 x 10° 2.32 x 10"^ 9.40 x 10"^ 

1 8.42 1.55 X 10^ 6.65 x 10"^ 2.19 x lO'^ 4.82 x 10^^ 
1 5.59 3.47 X 10^ 3.24 x 10"^ 3.38 x 10"^ 1.42 x lO""* 
5 9.50 7.90 x 10^ 2.32 x lO^^ 1.05 x lO'^ 4.20 x lO"*^ 

15 4.59 7.13 X lO^^ 1.43 x 10° 1.54 x lO'^ 2.61 x 10"^ 

1 8.39 3.15 x 10^ 3.18 x 10"^ 6.32 x 10"^ ~ 

1 5.23 5.64 X 10^ 5.02 x 10"^ 6.05 x lO"'' 3.19 x 10"^ 

5 8.66 1.37 x 10^ 7.35 x lO^^ 4.82 x 10"^ ~ 

30 9.15 4.43 x 10^ 2.30 x lO^^ 3.88 x 10"^ ~ 

30 6.62 9.75 x 10^ 1.42 x 10"^ 3.93 x 10"^ ~ 

60 9.52 6.31 x 10^ 1.63 x lO'^ 4.40 x 10"^ ~ 

60 7.03 9.44 x 10^ 1.45 x 10"^ 3.92 x 10"'' ~ 

600 11.04 9.68 x 10^ 6.84 x 10"^ 6.73 x lO''^ ~ 

3600 24.28 5.48 x 10^ 7.68 x 10"^ 5.86 x 10"^ ^ 

3600 16.47 9.99 x 10'^ 5.89 x 10"^ 4.89 x lO''^ ~ 

3600 13.00 3.40 x 10^ 2.40 x 10"^ 2.11 x 10"^ ~ 

20 1 10.03 3.61 X 10^ 2.94 x lO'^ 1.65 x 10"^ 1.95 x IQ-'^ 
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Table 2— Continued 



Mass (Mq)1 t (s)2 Period (s)^ t (s)"* l/r„„, (s-i)^ -1/r, (s-i)^ ^l/r,eafe (s-^)" 

1 6.77 1.77 X 10^ 9.49 x lO^^ 3.68 x 10"^ 1.43 x 10"'' 

5 11.14 2.91 x 10^ 4.18x10-3 7.44 x 10"^ 5.11 x 10"*^ 

25 5.86 7.29 x 10"! 1.38 x 10° 4.03 x 10"^ ^ 

1 10.70 2.89 X 10^ 3.79x10-2 3.27 x lO^^ ~o 

1 7.09 1.74 xlO^ 5.96x10-2 2.20 x IQ-^ ~0 

1 6.74 1.17 xlO^ 1.34x10-2 4.85 x 10-^ ~o 

5 11.30 3.97 xl02 3.91x10-3 1.39 x 10"^ ~0 

5 7.88 2.54 xl02 5.04 x 10"^ 1.10 x lO^^ ~o 

30 8.71 2.59 X 10^ 1.12 x 10"^ 7.33 x lO"" ~0 

60 9.18 1.01x10" 7.30x10-4 6.31 x 10"" ~0 

3600 17.25 8.35 x 10^ 6.13 x 10-^ 6.01 x 10"^ ~ 

30 5.35 2.28x10° 4.41 x 10-^ 3.31 x IQ-^ ~0 

1 7.51 2.44 xlQi 4.18x10-2 8.26 x lO"" ~0 

5 8.67 4.23 x 102 3.04 x 10-^ 6.76 x lO"" ~ 

30 9.79 3.84 x 10^ 7.83 x 10-" 5.23 x IQ-" ~ 

3600 16.13 5.74 x 10° 3.21 x 10-^ 3.20 x 10"^ ~ 

35 5.09 3.49 x 10^ 9.27 x IO-2 6.41 x IO-2 1.49 x 10-^^ 

4.82 5.28 X 10-^ 1.90 x 10° 6.93 x 10"^ 1.08 x 10-^ 

1 7.99 1.81 X 102 8.57 x 10-^ 3.05 x 10"^ -. 
1 7.60 1.35x102 1.06x10-2 3.16 x IQ-^ ~0 
1 5.85 5.22 X IQi 2.14x10-2 2.20 x IQ-^ ~0 
5 8.36 9.90x102 2.06 x 10-^ 1.05 x IQ-^ ~0 
5 6.44 1.09 xlO^ 2.14x10-3 1.22 x IQ-^ ~0 

30 9.19 5.03x103 7.07x10-* 5.09 x 10"" ~0 

60 9.70 6.22 x 103 5.61 x 10-" 4.00 x IQ-" ~ 

40 5.95 5.45 x 10* 6.82 x IO-2 4.98 x IO-2 1.83 x 10"^ 

5.47 5.43 X 10-^ 1.85 x 10° 7.20 x 10-3 1.55 x 10-3 

1 9.00 7.88 X 10* 1.60x10-2 3.31x10-3 ~0 
1 8.56 3.63 X 10* 3.03x10-2 2.81x10-3 ~0 
1 7.39 4.46 X 103 4.74 x 10-3 4.51 x 10-3 ~ 
1 6.36 5.41 X 10* 2.24 x IO-2 3.47 x 10-3 5.00 x 10-" 
1 5.95 5.62 X 10* 2.23 x IO-2 3.36 x 10-3 1.11 x 10-3 
5 9.45 6.06 X I02 2.55 x 10-3 8.97 x IQ-" - 

£ = 2 

11 1.36 4.55 X 10* 1.19 x 10-* 8.66 x IO-2 9.96 x 10-3 

1 4.47 9.37 x 102 6.07 x 10-3 5.00 x 10-3 ~ 
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Table 2 — Continued 



Mass (Mo)i t (s)2 Period (s)^ r (s)^ l/r„„, (s~i)5 -l/r, (s-i)^ -I/t,,,^ (s-i)^ 

1 2.63 1.02x102 1.57x10^2 5.87 x 10"^ 1.19x10"^ 

1 1.90 1.28 x 10^ 1.75 X 10"^ 7.99 x 10"^ 1.67 x 10-^ 

1 1.73 7.23 x 10^ 3.58 x 10"^ 4.37 x 10"^ 1.76 x 10"^ 

5 4.82 2.40 X 10^ 2.28 x lO'^ 1.86 x lO'^ ~ 

5 2.98 1.98 x 10^ 6.66 x lO^^ 1.59 x lO'^ 1.11 x 10"^ 

5 2.68 6.02 x 10^ 2.93 x 10-^ 2.74 x 10"^ 2.73 x 10"^ 

5 2.21 5.67 X 10^ 4.89 x 10"^ 2.40 x lO'^ 7.25 x lO""* 

5 1.98 1.41 x 10^ 1.15 X 10-2 1.31 x 10"^ 9.52 x lO^^ 

30 5.35 3.45 x 10^ 3.77 x lO'^ 8.73 x lO'^ ~ 

30 3.77 3.43 x 10^ 1.70 x 10"^ 1.41 x lO'^ ~ 

30 3.44 8.76 x 10^ 1.21 x 10"^ 7.12 x lO"" ~0 

30 3.07 6.34 x 10^ 2.98 x lO'^ 1.39 x 10"^ 1.02 x lO-^ 

30 2.55 2.50 x 10^ 5.46 x lO'^ 1.11 x lO'^ 3.54 x 10"^ 

30 2.29 2.94 x 10^ 6.24 x lO'^ 9.68 x 10"'' 1.87 x lO'^ 

60 5.74 6.82 x 10^ 2.08 x 10"^ 6.13 x lO"'' ~ 

60 4.10 2.49 x 10^ 1.30 x 10"^ 8.98 x lO"" ~ 

60 3.82 1.81 X 10^ 6.04 x lO'^ 5.00 x lO"" ~ 

60 2.78 1.21 X 10'"' 1.79 x 10"^ 7.59 x lO"'' 2.04 x 10~* 

60 2.50 6.38 x 10^ 1.79 x lO^^ 7.31 x lO"'' 8.98 x lO"'' 

13 2.68 2.77x10° 4.87 x lO'^ 1.26 x lO'^ 9.02 x 10^^ 

2.65 5.43 X 10"! 1.88 x 10° 3.80 x lO'^ 4.10 x 10"^ 

1 4.99 1.49 X 10^ 6.93 x lO'^ 2. 10 x lO'^ 1.89 x 10"^ 
1 3.44 3.46 x 10^ 3.24 x 10"^ 3.41 x lO'^ 6.10 x 10"^ 
5 5.64 7.10 X 10^ 2.41 x lO'^ 9.96 x lO"'' 2.04 x 10"^ 

15 2.85 6.14 X lO^^ 1.64 x 10° 1.25 x 10"^ 7.38 x lO^^ 

1.79 1.01 X 10^ 1.76 X 10^1 6.39 x 10"^ 1.35 x 10"^ 

1.72 3.76 X 10° 3.35 x lO'^ 3.67 x 10"^ 3.26 x 10"^ 

1 6.16 5.42 X 10^ 1.85 x lO'^ 3.28 x 10"^ ~ 
1 3.24 1.89 x 10^ 6.47 x 10"^ 4.8O x lO"'' 1.14 x 10"^ 
5 6.34 1.75 X 10^ 5.76 x lO'^ 2.55 x lO''^ ~ 

60 4.34 1.10 x 10^ 1.24 x 10"^ 3 28 x 10"^ ~ 

3600 14.17 6.48 x 10'* 7.16 x 10"^ 5.62 x lO''^ ~ 

3600 10.32 1.58 x 10^ 2.71 x 10"^ 2.08 x 10"^ ~ 

3600 9.11 1.73 xlO^ 4.80 x 10"^ 4.22 x lO''^ ~0 

20 1 6.04 3.50 X 10* 3.02 x 10"^ 1.65 x 10"^ 9.49 x lO''^ 

1 4.19 3.43 x 10^ 7.12 x lO'^ 4.15 x 10"^ 4.97 x 10"^ 

25 2.38 2.85x10° 4.67x10"* 7.70x10"'' 1.16x10"* 

1 6.43 3.07x10* 3.59 X 10"2 3.39 x 10"^ ~0 
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Table 2— Continued 



Mass (Mq)1 t (s)2 Period (s)^ t (s)"* l/r„„, (s-i)^ -1/r, (s-i)^ ^l/r,eafe (s-^)" 

1 4.43 1.52x101 6.78x10-2 1.87 x lO^^ ~0 

1 4.24 2.20 x 10^ 9.57 x lO'^ 5.02 x 10"^ ^ Q 

5 6.58 1.21 X 10^ 2.49 X 10~3 1.67x10^3 ^Q 

5 4.91 2.32 X 10^ 5.39x10^3 1.07x10^3 ~0 

30 5.41 1.56 xlO^ 1.17x10-3 5.28 x lO"" ~0 

60 5.70 8.32 x 10^ 7.60 x 10"^ 6.40 x lO"" ~ 

3600 4.59 4.81 X 10^ 3.23x10"^ 1.15 x 10"^ ~0 

30 3.20 2.27x10° 4.41x10^1 1.12x10^3 ~o 

1.68 1.49x10° 7.94x10-1 2.55x10^3 1.21x10^1 

1 4.52 2.47 xlO^ 4.12 x 10"^ 6.75 x 10"^ ~0 
5 5.23 4.11 X 10^ 3.01x10-3 5.74 x 10"^ ~0 

30 5.91 3.37 X 10^ 7.55 x 10"^ 4.58 x 10"* ~ 

35 2.98 4.93 x lO^i 2.03 x 10" 5.28 x 10"^ 6.57 x lO"'! 

1 4.77 2.31 X 10^ 7.53x10-3 3.20x10^3 ^o 

1 4.56 1.34 xlO^ 1.06x10-2 3.15 x 10-^ ^o 

1 3.63 5.13 X 10^ 2.17x10-2 2.23 x 10-^ ~0 

1 3.41 2.01x103 4.80x10-3 4.30x10-3 ~0 

5 5.00 8.02 X 10^ 2.24 x 10-3 9.89 x 10-" ~ 

5 3.99 1.45x103 1.96x10-3 1.27x10-3 ~0 

30 5.49 4.08 x 103 7.37 x 10"" 4.92 x lO"" ~ 

60 5.79 4.78 x 103 5.93 x lO"* 3.84 x lO-'' ~ 

40 3.44 5.12 x IQ-^ 1.96 x 10° 6.00 x 10-3 7.37 x lO-"* 

2.04 1.14 xlO^ 2.73x10-1 1.27x10-3 1.83 x 10-^ 

1 5.37 7.42x101 1.68 x 10-^ 3.32x10-3 ~0 
1 5.14 3.53x101 3.11x10-2 2.85x10-3 ~0 
1 3.99 6.50 X lOi 1.94 x 10-2 3.74 x 10-3 2.70 x 10-" 
1 3.75 4.21 X lOi 2.75 x 10-2 2.95 x 10-3 8.16 x 10"" 
5 5.68 5.29x102 2.75x10-3 8.57x10-" ~0 

11 0.99 1.10 xlOi 1.63x10-1 5.81x10-2 1.35x10-2 

0.88 3.74 x 10° 4.71 x lO-i 1.03 x 10-2 1.94 x lO-i 

1 1.95 1.12 X 102 1.42 X 10-2 5.20 x 10-3 2.74 x 10-^ 
1 1.33 2.69 x IQi 4.67 x IO-2 3.09 x 10-3 6.43 x 10-3 
5 3.49 1.49 x 103 2.34 x 10-3 1.66 x 10-3 ~ 
5 2.21 2.90x102 5.19x10-3 1.73x10-3 1.77 x 10-^ 
5 2.04 7.05 x 102 3.85 x 10-3 2.43 x 10-3 6.93 x 10-^ 
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Table 2 — Continued 



Mass (Mo)i t (s)2 Period (s)^ r (s)^ l/r„„, {s-^f -l/r, (s^^f -I/t,,,^ (s-y 

5 1.72 3.69 X 10^ 2.98 x 10"^ 2.63 x 10"^ 8.16 x 10"^ 

5 1.53 9.65 x 10^ 1.46 x 10"^ 1.06 x 10"^ 3.19 x IQ-^ 

30 3.86 3.89 x 10^ 3.45 x IQ-^ 8.84 x lO""^ ~ 

30 2.54 8.66 x 10^ 1.22 x 10"^ 6.88 x lO"" ~ 

30 2.34 3.47 x 10^ 4.18 x 10"^ 1.30 x 10"^ ^ 

30 1.98 3.78 x 10^ 3.92 x IQ-^ 1.22 x 10"^ 4.68 x 10"^ 

30 1.79 1.43 X 10^ 8.31 x 10"^ 8.23 x lO^'' 4.93 x lO"-* 

60 4.15 8.29 X 10^ 1.86 x IQ-^ 6.54 x lO"" ~ 

60 2.83 1.53 x 10^ 6.97 x lO'^ 4.27 x 10"^ ~ 

60 2.56 5.36 x lO" 9.63 x IQ-^ 9.45 x 10"'' ~ 

60 2.16 2.59 x 10^ 1.23 x IQ-^ 8.23 x 10"" 2.29 x IQ-^ 

60 1.95 7.26 X 10^ 2.24 x IQ-^ 6.45 x 10"'' 2.18 x lO"* 

13 2.02 9.43 x IQ-i 1.14 x 10° 8.20 x 10"^ 3.31 x IQ-^ 

2.01 9.39 X 10-1 1.15 x 10° 8.29 x 10"^ 3.63 x 10"^ 

2.01 9.39 X 10-1 1.15 x 10° 8.29 x 10"^ 3.63 x lO"*^ 

1 3.66 1.41 X IQi 7.29 x 10"^ 1.99 x lO'^ 3.30 x 10"^ 
1 2.65 3.66 X IQi 3.09 x 10"^ 3.53 x 10"^ 7.66 x 10"^ 
5 4.14 6.59 x 10^ 2.47 x 10"^ 9.52 x 10"'' 3.54 x lO"'^ 

15 2.21 5.65 X lO^i 1.77 x 10° 1.06 x 10"^ 1.26 x lO^^ 

1.42 3.08 X 10^ 1.02 x 10^^ 6.84 x lO'^ 1.45 x lO'^ 

1.38 4.50x10° 2.57x10-1 3.07 x lO'^ 4.45x10^3 

1 2.49 1.31 X IQi 7.92 x 10^^ 3 33 ^ IQ-^ 2.56 x IQ-^ 
1 1.67 2.17x101 5.90x10-2 1.34x10-3 1.16 x 10"^ 
5 2.76 1.89 X 10^ 7.16 x IQ-^ 2.20 x 10"'' 1.64 x lO^^ 

30 3.06 1.05 x 10^ 1.22 x IQ-^ 2.67 x 10"'' ~ 

60 3.26 1.07 xlO^ 1.20 x 10"^ 2.67x10-4 ~0 

3600 10.18 8.41 X lO"! 6.51 x 10^^ 5.32 x IQ-'^ ~ 

3600 8.42 2.73 x 10^ 8.11 x 10^^ 4.45 x 10"^ ~ 

3600 6.95 5.68 x 10^ 5.03 x 10"^ 4.85 x IQ-'^ ~ 

20 1 4.58 6.79x101 1.79 x 10"^ 3.14 x lO'^ 7.57 x lO"'^ 

1 4.51 6.96 x lOi 1.76 x 10"^ 3.25 x 10"^ 8.71 x 10"^ 

1 3.35 2.38 X 10^ 8.13 x 10"^ 3.91 x lO'^ 5.18 x 10^^ 

25 2.00 2.80 x 10° 3.82 x lO^i 3.23 x 10"'' 2.42 x 10"^ 

1.00 1.17 X 10^ 1.99 X 10-2 9.07 x lO'^ 2.28 x lO'^ 

1 3.45 1.21 X IQi 8.42 x lO'^ 1.48 x 10"^ - 
1 2.20 2.17 xl02 2.06x10-2 2.23 x lO'^ 1.37 x 10"^ 
1 2.17 3.40 X 10^ 1.74 x lO'^ 2.94 x lO'^ 1.15 x 10"^ 
5 3.81 1.92 X 10^ 6.24x10-3 1.03x10-3 ~0 
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Table 2 — Continued 



Mass (Mo)i t (s)2 Period (s)^ r (s)^ 1/t„„, {s-^f -l/r, {s'^ -l/ri.ak (s-^Y 

30 4.20 2.15 X 10^ 1.24 x lO^^ 7.70 x lO""^ - 

60 4.41 9.52 X 10^ 7.75 x lO"'' 6.70 x lO"" ~ 

600 3.18 2.42 X 10^ 4.36 x 10"^ 2.25 x 10"^ ~0 

3600 3.86 4.35 x 10'' 3.23 x 10"^ 9.29 x 10"^ ~ 

30 2.36 2.30 x 10° 4.35 x IQ-i 6.71 x 10"^ ~ 

1.38 1.10x10° 9.31x10-1 1.52x10^3 2.41 x 10"^ 

1 3.34 2.44 xlO^ 4.16 x 10"^ 5.30 x lO"" ~0 
1 1.66 6.45 X 10^ 4.43 x IQ-^ 4.29 x 10^^ 2.83 x 10"^ 
5 3.87 3.90 X 10^ 3.01 x 10"^ 4.51 x 10"^ ~ 

30 4.38 2.59 x 10^ 7.54 x IQ-^ 3.68 x 10"" ~ 

35 2.30 4.93 x 10"^ 2.03 x 10° 2.93 x 10"^ 8.37 x lO"'^ 

1 3.55 6.73 X 10^ 4.07 x 10"^ 3.92x10^3 ^0 

1 3.41 1.11 X 10^ 1.20x10-2 2.98x10-3 ^0 

1 2.79 4.35 X 10^ 2.49 x 10"^ 1.93 x lO^^ ~o 

1 2.65 2.13 X 10" 4.25x10-3 4.20 x 10-^ ^Q 

5 3.71 7.14 xlO^ 2.34x10-3 9.42x10-" ~0 

5 3.07 1.42x103 1.94x10-3 1.23x10-3 ~0 

30 4.06 6.22 x 103 6.79 x IQ-" 5.18 x 10"" ~ 

60 4.28 6.63 x 103 5.36 x IQ-" 3.85 x 10-" ~ 

3600 3.31 1.63x10'"^ 1.14x10-^ 5.26 x 10-^ ~0 

40 1.74 4.63 x 10° 2.59 x IQ-^ 3.73 x 10-" 4.29 x IQ-^ 

0.98 2.71 X 10^ 2.07x10-2 1.51 x 10-^ 1.89x10-3 

1 4.00 1.24 X 10^ 1.20 X IQ-^ 3.90 x 10-3 ~ 
1 3.86 2.85 X 10^ 3.77 x IQ-^ 2.67 x 10-3 - 
1 3.55 2.18 X 103 5.41 x 10-3 4.95 x 10-3 5.42 x IQ-'^ 
1 3.09 4.21 X 10^ 2.70 x IQ-^ 3.22 x 10-3 6.66 x IQ-^^ 
1 2.94 4.82 X 10^ 2.41 x 10-^ 3.23 x 10-3 1.16 x IQ-" 
5 4.25 4.14 X 10^ 3.18 x 10-3 7.66 x 10-" - 
5 3.20 3.03 X 103 1.66 x 10-3 1.22 x 10-3 1.08 x IQ-" 

30 4.65 1.90 X 10" 4.75 x 10-" 4.23 x 10-" - 

^The initial mass of the model. 

^The time before collapse in seconds. 

3 The period of the particular growing mode. 

"The growth timescale given by 1/r = l/r„„c + 1/t^ + 1/T/eafc 

^Inverse of the growth timescale from nuclear burning. 



-20- 

^Negative and inverse of the decay timescale due to v emission. 

^Negative and inverse of the decay timescale due to acoustic losses to the envelope. 

®Each row is associated with a specific i given by the header for each section. 
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Fig. 1. — Density (top panel) and electron fraction, Ye, (bottom panel) versuses interior mass, M, for the 11 
Mq (solid and dark blue), 13Mq (dotted and light blue), 15 M0 (dashed and green), 20 M0 (dot-dash and 
olive green), 30 Mq (dot-dot-dot-dash and red), 35 M0 (long dashed and orange), and 40 M© (solid orange) 
models, all at 5 s before the onset of collapse. Note in the top panel that the lower mass models have higher 
central densities, but that the higher mass models have shallower density gradients. Included on the Ye plot 
are points indicating significant features in the star, which greatly affect the modes. Circles represent the 
farthest extent of the fossil Fe core, which was formed during convective core Si burning. Stars represent 
the base of the Si burning shell at 5 seconds before the start of collapse (i.e., the extent of the Fe core, 
including the fossil Fe and the most recent ashes from Si shell burning), which generally aren't convective 
until collapse. The triangles represent the base of the convective O burning shell 
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Fig. 2. — Density (solid green) and Ye (solid red) versuses interior mass in comparison to the propagation 
diagram for the 13 Mq model at 5 s before the onset of collapse. Usually, the square of the Brunt- Vaisala (iV^) 
and Lamb (Lf) frequencies are plotted. In this case, an equivalent period of oscillation is displayed instead, 
defined by Psrunt = 27r/-\/(|A^^|) and Phamb = 27r/-\/(|L^|). The circles, stars, and triangles represent the 
same structural boundaries as in Fig. 1. The vertical spikes are from discontinuous density and Y^ changes 
in the profiles and are real. The bottom panel provides the equivalent information for the 30 Mq model. 
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Fig. 3. — A time series of the propagation digram for the 13 Mq model from 600 s to s before the start of 
collapse. Psrunt and Phamb are defined as in Fig. 2. Notice that the evolution aburptly changes for t = s. 
This is due to the dramatic increase in nuclear rates during implosive burning, which causes new regions to 
become convective in the mixing-length calculation. 



24 



1 - 



"D 
O 

*^ 
CD 

Q_ 



10 



. M = 15M0 


> 


1 1 1 1 

N 




t = 5 s: 

t = ^ : 










\ 






Ix 




\ 


" 


1 


t 


\ 


" 




}\k, 




^^ 






^^^* 




1 


^^^^ 






^^^ 


1 


1 . . . . 


1 





0.5 



1.0 
M (Mo) 



1.5 



Fig. 4. — The propagation diagram for ^ = 1 and the 13 Mq model at 5 s before the onset of collapse with 
Psrunt (solid) and PLamb[(- = 1) (dashed), as defined in Fig. 2. The horizontal Hnes correspond to each 
eigenmode and period and are quantized. The shading in the lines is a measure of the magnitudes of the 
square of the amplitude, |^p. Within a single line brighter portions indicate a higher amplitude while the 
darker color signifies a lower amplitude. The blues indicate modes which are stable (r < 0), and the reds 
correspond to the modes which are unstable (r > 0). 
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Fig. 5. — This figure is similar to Fig. 4, but for the 30 M© model. 
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Fig. 6. — The dimcnsionlcss cigenfunctions of stable g-modc oscillations with periods of 10.51, 9.28, 8.17, 
and 7.07 s, for ^ = 1, 13 M©, and 5 s before the start of collapse. Each mode is distinguished by its period. 
The top panel plots S,r / R, where ^^ is the radial component of perturbation and R is the calculation domain 
size. The amplitude is arbitrarily scaled to 1.0 at its peak value. The bottom panel plots the corresponding 
horizontal amplitude, ^h/R- 
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Fig. 7. — Plotted are the same attributes depicted in Fig. 6, but for a representative unstable mode with 
period 9.50 s. 
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Fig. 8. — The Lagrangian displacement squared, |^p (browns and oranges), the local deposition of energy 
by nuclear processes, e„„c (solid and black), and the local losses of energy due to neutrinos, €i, (dashed), 
versus interior mass (Mq) for the 13 Mq model, ^ = 1, and 5 s before the onset of collapse. The top panel 
illustrates this comparison for two representative stable modes (periods — 10.51 and 9.28 s). while the bottom 
panel plots an unstable mode (period = 9.5 s). Note that the inner-core g-modes have their largest relative 
amplitudes where neutrino losses dominate, but that the outer-core g-modes have their largest relative 
amplitudes trapped in the recent Si burning ashes as well as the Si burning layer giving rise to instability. 
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Fig. 9. — Same as Fig. 8, but for the AI = 30 M© progenitor, stable modes with periods of 10.48 s and 9.19 
s, and an unstable mode with period = 8.67 s. 
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Fig. 10. — The growth timcscales, r, in seconds, for the unstable modes versus their periods, in seconds, for 
the 13 Mq model. The plots are for £ = 1 (circle), 2 (triangle), and 3 (square) and for times of 60 (blue), 
30, 25, 20, 15, 10, 5, 2, 1, and (red) seconds before the onset of collapse. The arrow indicates the evolution 
of the period and growth timcscales for particular modes as the models approach collapse. 
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Fig. 11. — Same as Fig. 10, but for the AI = 30 Mq progenitor. Note the clear, unintcrupted evolution of 
specific g-modes up to the onset of collapse. 
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Fig. 12. — Same as Fig. 10 but for the i\/ = 11 M0 progenitor. While plenty of growing modes exist, their 
timcscales are far too long for any significant growth to occur. 



